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Abstract 

We analyse the experimental time series of internal pressure in a four cylinder spark ignition 
engine. In our experiment, performed for different spark advance angles, apart from usual cyclic 
changes of engine pressure we observed oscillations. These oscillations are with longer time scales 
ranging from one to several hundred engine cycles depending on engine working conditions. Basing 
on the pressure time dependence we have calculated the heat released per cycle. Using the time 
series of heat release to calculate the correlation coarse-grained entropy we estimated the noise 
level for internal combustion process. Our results show that for a smaller spark advance angle the 
system is more deterministic. 
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A combustion process in spark ignition engines is known as nonlinear and noisy 
one. Combustion instabilities are occurring as a cycle-to-cycle variations of in- 
ternal cylinder pressure effecting directly on the power output. Examination 
of these variations can lead to better understanding of their sources and help 
in their eliminations in future. Improving engine efficiency requires achiev- 
ing better combustion conditions without introducing additional disturbances. 
In the present paper we analyse the dynamics and estimate the noise level in 
combustion process basing on experimental time series of internal pressure and 
calculated from them heat release. In the following analysis we apply the non- 
linear multidimensional methods which can distinguish random variations from 
a deterministic behaviour. 



I. INTRODUCTION 



Combustion in four stroke spark ignition (SI) engines is a complex cyclic process con- 
sisted of air intake, fuel injection, compression, combustion, expansion and finally gas ex- 
haust phases (Fig. |IJ) where burned fuel power is transmitted through the piston to the 
crankshaft. In early beginning of SI engine development there were observed instabilities 
of combustion 1 . These instabilities are causing fluctuations of the power output making it 
difficult to control. The problems of their sources identification and their elimination have 
became the main issues in SI engines technologies engineering and they have not been solved 
up to present time^. Among the the main factors of instabilities classified by Heywood^ are 
aerodynamics in the cylinder during combustion, amounts of fuel, air and recycled exhaust 
gas supplied to the cylinder and a local mixture composition near the spark plug. 

Recently, Daw et al&L and Wendeker et a/.— have done the nonlinear analysis of such pro- 
cess. Changing an advance spark angle they observed the considerable increase of pressure 
fluctuations level^ claiming that it is due to nonlinear dynamics of the process. In the other 
work^ Wendeker and coworkers have proposed intermittency mechanism to explain the rout 
to eventually chaotic combustion. 

Prompted by these findings we decided to analyse the the correlation entropy of the 
combustion process in different working conditions of the engine. With a help of entropy 
produced by the dynamical system we can quantify the level of measurement or dynami- 
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FIG. 1: Schematic picture of a combustion cycle in a four stroke spark ignition engine. 

cal noise^. In the present paper we shall start our analysis from examining experimental 
pressure time series. 

It should be noted that pressure is the best known quantity to analyse engine dynamics. 
Cylinder pressure together with volume data can be used to obtain indicated mean effective 
pressure (IMEP), calculate the engine torque, indicated efficiency and also burn rate, bulk 
temperature and heat release. Moreover, statistical analysis of the pressure data can also 
provide information about combustion process stability. 

However, in practice, It is not easy to perform direct measurement of pressure^, as one 
needs a good sensor persistent to hight temperatures to be placed inside the engine cylinder. 
Therefore to obtain information about pressure some researchers developed alternative non- 
direct measurements procedures^. 

In our case we have been dealing with novel pressure fibre optical sensors^. Due to 
applying them noise from measurement is very low, comparing to traditional piezo-electric 
ones^. This enabled us to examine the dynamics more effectively than it was possible in 
earlier investigations. 

The present paper is divided into 6 sections. After present introduction (Sec. 1), we 
will provide the description of our experimental standing and measurement procedure in 
Sec. 2. There we present some examples of cycle-to-cycle variations in pressure inside one 
of cylinder. In Sec 3. we examine the pressure with more detail. We also perform spatio- 
temporal analysis comparing the fluctuations of pressure in succeeding cycles for different 
advance angles. In Sec. 4 we calculate the heat release per cycle. Finally in Sec. 5 we 
analyse its time dependence and show our main result i.e. level of noise. We end up with 
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FIG. 2: Experimental stand. 



conclusions and last remarks (Sec. 6). 



II. EXPERIMENTAL FACILITIES AND MEASUREMENTS OF INTERNAL 
PRESSURE 

In our experimental stand (Fig. EJ) pressure was measured directly inside cylinder by the 
use of the optical fibre sensorJA. Such equipment provides one of the most direct measures of 
combustion quality in an internal combustion engine. Internal pressure data were obtained 
from Engine Laboratory of Technical University of Lublin, where we conducted series of 
tests. 

The pressure traces were generated on a 1998 cm 3 Holden 2.0 MPFI engine at 1000 RPM. 
The data was captured by use of NuDAC-TK v. 2.0 data acquisition and data processing 
program 15 . The original files contained cylinder pressure at crank angles 0-720 degrees. Each 
of three large file (about 990 MB per each) contained above 10000 combustion cycles. Data 
was taken at different spark timings (spark advance angles): 5,15,30 degrees before top dead 
center (BTDC). The engine speed, air/fuel ratio, and throttle setting were all held constant 
throughout the data collection period. Intake air pressure, in inlet pipe, was value 40 kPa. 
Torque for each of three spark timings were adequately: 21, 28 and 30 Nm. 

To perform signal analysis we needed large enough data. In this aim we measured 10000 
cycles for each of three spark advance angles Aa z . The results for first 1000 cycles are shown 
in Figs. 3a-c. Note that depending om an advance angle we have more or less broadened 
region of pressure fluctuations. The full line shows the pressure averaged over the first 1000 
cycles. It is increasing with growing Aa z and reaches its highest value for Aa z = 30°. 

Note also , in our for four stroke engine the combustion period (Figs. 1,3) corresponds 
exactly to the double period of the crankshaft revolution synchronised with a single spark 
ignition. Every combustion cycle starts with initial conditions given by a mixture of air and 
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TABLE I: Definitions of variables and symbols used in the paper. 



position of the piston 


h 


internal cylinder pressure 


P 


actual cylinder volume 


V{a) 


Heaviside step function 


G(z) 


heat released 


Q 


heat released in particular cycle i 


Qi 


heat released vector in embedding space 


Q 


spark advance angle 




Aa z embedding time delay in cycles 


m 


cycle number 


i,j 


embedding dimension 


n 


number of considered points in time series 


N 


loading torque 


F 


crank angle 


a G [0, 720°] 


threshold 


£ 


correlation integral 


c n 


coarse-grained correlation integral 


C n {e) 


correlation entropy 


K 2 


coarse-grained correlation entropy 


K2(e) 


calculated from time series 




coarse-grained entropy 


Knoisy 


correlation dimension 


D 2 


Noise-to-Signal ratio 


NTS 


standard deviation of data 


&DATA 


error function 


Erf(z) 


fitting parameters 


X , a,b 


standard deviation of noise 


a 



5 



0- 0,8 
^0,6 
3 0,4 
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FIG. 3: Internal pressure of 1000 combustion cycles against a crank angle for a spark advance 
angle Aa z = 5, 15 and 30 degrees for Figs. 3a, b and c, respectively. Full lines correspond to the 
average angular pressure. 

fuel. All of succeeding cycles are separated by gas exhaust and intake stroke phases. That 
process is in general nonlinear and can be also mediated by stochastic disturbances coming 
from e.g. non homogenous spatial distribution of fuel/air ratio. After combustion exhaust 
gases are mixed with fresh portions of fuel and air. 

Therefore the residual cylinder gases after each combustion cycle influence the process in 
a succeeding cycle leading to different initial conditions of air, fuel and residual gas mixture 
contents. 
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TABLE I: continuation 



cylinder diameter 


D 


= 86 mm 


crank radius 


r 


= 43 mm 


connecting-rod length 


I = 


= 143 mm 


heating value of the fuel 


W u = 


43000 kJ/kg 


compression ratio 




e = 8.8 


Poisson constant 


K - 


= 5£ « 1.4 


mass burned in cycle i 




Mi 


autocorrelation function 




AC(j) 


output torque 




S 



ILL ANALYSIS OF PRESSURE 



During the combustion process the internal volume of engine cylinder is driven kinemat- 
ically by the piston. In a consequence of above the internal it changes also periodically as a 
function of crank angle a and satisfying the relation 



D 2 D 2 1 
V(a) = K—^-h + TT—2r -. 



(1) 



where the piston position h 




(2) 

and constants r, I, D as well defined in Tab. |l] 

In some sense the combustion initiated by ignition in each engine cycle is an indepen- 
dent combustion event. Such events are separated by the processes of exhaust and intake 
dependent on the amount of combusting fuel mass and quality of newly prepared fuel-air 
mixture. To illustrate this effect we are showing in Figs 4a-c spatio-temporal plots corre- 
sponding to first 1000 cycles of our pressure time series for different advance angle Aa z = 5, 
15, 30 degrees, respectively. Each colour in Figs. 4a-c correspond to one of four interval of 
[P min , PmaxY white [P min , P^, green [Pi,P 2 ], red [P 2 ,-Ps] and blue [P^P m ax\- One can easily 
see that the pressure signal, especially in Fig. 4c seems to change in some regular manner of 
a time scale of about 100 cycles. Similar feature is also visible in Fig. 4a while it is difficult 




Number of cvcle 



(c) 



FIG. 4: Four colour Spatio-temporal corresponding combustion process parameters as in Fig. 
3a-c Respectively. Each colour correspond to one of four interval of [P m i n , Pmax\- white [P m i n , P\], 
green [Pi,P 2 ], red [P 2 ,P 3 } and blue [P 3 ,P max } (Pmin = -0.2MPa, Pi = O.lMPa, P 2 = 0.2MPa, 
P 3 = 0.6MPa, P max = 2.0MPa). 

to find such regularity in Fig. 4b. However after more careful examination one can identify 
such a time scale consisting of about 450 cycles. 

The broad angular interval of fluctuations visible in Fig. 3b. has its consequences in 
irregularly border between red and blue colours (Fig. 4b). 



IV. VARIATIONS OF HEAT RELEASE 

To capture the cycle-to-cycle changes in combustion process we decided to calculate of 
heat release in a sequence of combustion cycles. This quantity is more convenient to examine 
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FIG. 5: Heat release per cycle versus sequential cycles. 



stability of combustion process because enables us to concentrate on it. In contrast to that 
internal pressure was is effected by combustion and cyclic compression. Heat release has 
also practical meaning as it is proportional to burned fuel mass. 

For an adiabatic process the heat released from chemical reactions during combustion in 
the engine is given by an differential equation coming from the first law of thermodynamics 
with respect to a crank angle a: 



dQ k dV 1 Tr dp 

— = p 1 V — . 

da k — 1 do: k — 1 da 



(3) 



Using it together with the pressure time series and parametric change of cylinder volume 
V Eqs. 1-2 we have done calculation of heat released in succeeding cycles Qi. 
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It is closely related to burned fuel mass in one cycle 



Mi = Qi/W, 



(4) 



where W u is heating value of the fuel, listed in Tab. 1. I should be underlined that the last 
equation (Eq. 4) neglects the effects of heat exchange between the cylinder chamber and its 
walls. This is in spirit of an adiabatic process assumption (Eq. 3). Of course the consumed 
mass will be larger because it also depends on the quality of mixture and combustion process. 

The calculated heat release Q i: for first 5000 cycles is plotted against cycles % in Fig. 5a-c 
for Aa z = 5, 15, 30 degrees, respectively. 

Note that in all cases there is some modulation ranges from one to few hundred cycles. 
Interestingly, for a small advance angles 5a z = 5° or 15° this modulation evolute indicating 
that the system can have quasiperiodic or chaotic nature. Note also for the first 1000 
variations for Aa z = 5° (Fig. 5a) resembles those for Aa z = 30° (Fig. 5c) while for 
Aa z = 15° the long time scale modulation is different. This is consistent with Fig. 4a-c. 
Generally, for Aa z = 30° the oscillations a of higher frequency and more regular. The high 
values of Qi in the middle part of Fig. 5c (i G [1300, 2400]) are connected with measurment 
instabilities. In that case however the average value of heat release < Qi > (< Qi >= 270J 
- Fig5a < Qi >= 237 J - Fig5b < Q { >= 18 1J - Fig5c) is the smallest indicating the lowest 
burning rate of fuel. In spite of that the output torque, for the same speed of a crankshaft, 
was relatively larger S = 30Nm (in the case Aa z = 30°) comparing to other levels 21 and 
28Nm for Aa z = 5° and 15°, respectively. Obviously, there are better combustion conditions 
in the last case. 

We have also calculated autocorrelation function from the whole 10000 cycles signal via 



with appropriate normalisation to one. The results for all three advance angles Aa z are 
depicted in Fig. 6. One can see that the decay of AC(j) amplitude with growing j is 
comparable but frequency of modulation is different for all these cases. Clearly, this is 
higher for larger Aa z . 



ACU) = "EQ(i)Q(i+j) 



(5) 
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FIG. 6: Correlation function calculated from heat release time series for different advance angle 
Aa z . 

V. ESTIMATION OF NOISE LEVEL FROM HEAT RELEASE SERIES 

In this section we shall examine the level of noise in heat release time series. In this aim 
we use nonlinear embedding space approach^. 

In the n dimensional embedding space the state is represented by a vector 

Qi+(n-l)m}, (6) 

where m denotes the embedding delay in terms of cycles. The correlation integral calculated 
in the embedding space can be defined a a 17 i 18 

i N N 

^) = ^EE e ( £ -IIQ*-Qili), (7) 

i j^i 

where N is the number of considered points corresponding to pressure peaks in cycles and 
is the Heaviside step function. For simplicity we use maximum norm. The correlation inte- 
gral C n (e) is related to the correlation entropy ^(e) and the system correlation dimension 
D2 by the following formul a 1 ^ 1 ^ 

lim C n (e) =D 2 \ne- nmK 2 (e). (8) 

n— >oo 

The coarse-grained correlation entropy can be now be calculated as 
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FIG. 7: Noise to signal ratio NTS versus a spark advance angle Aa z 



Ko(e) = lim In 



C n (e) 



d\nC n (e) 



(9) 



n ^oo c n+1 {e) dn 
In such a case the correlation entropy is defined in the limit of a small threshold e. 

In presence of noise described by the standard deviation o of Qj time series, the observed 
coarse-grained entropy K noisi j 10 ^ A can be written as 

1 



K n0 i S y\£) 



x 1 



e 2 /3 + 2a 2 - e/y/3 



(10) 



Function g(z), present in the above formula, reads 



2 ze 



(11) 



7T Erf (z) 

where Erf(.) is the Error Function. The parameters \i a ) b as well as a are unconstrained. 
They should be fitted in Eq. fllOl) to mimic the observed noisy entropy calculated from 
available data. 

After application the above method to the heat release times series we estimated noise 
calculating Noise to Signal ratio (NTS): 



NTS 



a 



O'DATA 



(12) 
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where <Jdata is the standard deviation of data. 

Figure 7 shows NTS versus a spark advance angle Aa z . It appeared that for any of 
examined cases the noise level is not high. Starting from a small spark advance angle 
Aa z = 5° the level of noise is the smallest NTS~ 13%. It grows to about 17% for Aa z = 15° 
and 25.5% for Aa z = 30°, respectively. In fact easy to note that all points (in Fig. 7) lay 
on a line, and ratio may be discribed as increasing linearly with Aa z . 

In all considered cases we have to do with a some large scale signal modulation ranging 
to a few hundered engine cycles and fluctuations of a few tens or more as well as very fast 
ones (Figs. 4 and 5). Using the same method we have examined the natue of very fast 
fluctuations and found that they are purely stochastic. Fluctuation of large time scales may 
involve in the estimation of noise level to make them underestimated by a few percents. 
This is because the correlated noise can occur and such a noise is irrelevant to our method. 
We think that in all three cases the calculated noise level are biased similary, because of 
the same experimental conditions, so the linear increasing behaviour of NTS versus advance 
angle Aa z is preserved. 

VI. REMARKS AND CONCLUSIONS 

In this paper we analysed instabilities of combustion process We started from analysing 
pressure time series. Using spatio-temporal methods we established that there is a long 
time scale in fluctuation of our experimental data. To examine this phenomenon in detail 
we calculated heat release and we performed the noise level estimation using nonlinear 
multidimensional methods. Our results clearly indicate that the noise in the time series is 
the highest for the largest advance angle. In case of Aa z = 30° we have got the signal with 
characteristic 100 cycles periodicity. 

Heat released in cycle as a practical parameter closely related to burned fuel mass and 
enables to follow the stability of combustion process better. Our noise estimation basing on 
heat release time series is more credible than the analysis pressure histories itself, as pressure 
is effected by volume cyclic compression and expansion phases. 

The method using correlation entropy which we applied here differs from the symbolic 
treatments'*^ used also for exploration of the engine dynamics by Daw et a& In their paper 
the signal was digitised and basing on the probabilities of 1 sequences probabilities the 
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information Shannon entropy was estimated. In our paper we estimate the NTS ratio by 
fitting coarse-grained entropy obtained from experimental time series to a general formula of 
correlation entropy evaluated in presence of noise. In our case the entropy has its dynamical 
meaning as a measure uncertainty of system state. 

One should also note that our present examination was limited to only one crank rota- 
tional speed 1000 RPM. Using the above procedure we now preparing systematic analysis 
with the other speeds for a future report. 
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